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ABSTRACT 


The orbiting Burst and Transient Source Experiment 
(BATSE) will locate gamma ray burst sources by analysis of 
the relative numbers of photons coming directly from a 
source and entering its prescribed array of detectors. In 

nyrlov fn af'f'iiral-Al v lnfat-fl hiiKst fioi l rrfts ih is thus 
v ^ — ~ ^ — A — — 

necessary to identify and correct for any counts contributed 

by events other than direct entry by a mainstream photon. 

The present work describes an effort to estimate the photon 

numbers which might be scattered into the BATSE detectors 

from interactions with the Earth's atmosphere. A model has 

been developed which yielded analytical expressions for 

single-scatter photon contributions in terms of source and 

r' -4 4- /-» 1 1 i 4- 1 r\or»4* i rvne 
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INTRODUCTION 


One of four experiment packages aboard the Gamma Kay 
Observatory (GRO) is the Burst Transient Source Experiment 
(BATSE). Utilizing an array of eight Nal detectors, BATSE 
will make angular location determinations of the gamma ray 
burst sources by means of analyzing the relative count 
responses of its detectors. Figure 1 depicts the GRO 
configuration 1 . If the planes of the detector faces were 
extended so as to form a closed surface, that surface would 
be the regular octohedron also depicted in Figure 1. 

The assignment of a numbering scheme to the detectors 
as indicated in Figure 2 caij be used to determine the 
angular coordinates, (9 , $ ), of a burst as follows. 

Suppose that, for a given burst, detectors 1, 2, and 3 are 
those presumed to optimally locate the burst source. Having 
recorded respective count rates N^, N 2 , and N 3 allows the 
(uncorrected for error sources, ^or " zeroth-order" ) angular 
coordinates to be computed from'*' 


V2 


tang* = [N x 2 + N 2 2 + N 3 2 + 2N 3 (N 2 - Nj)] , (1-a) 


tan<f> = (N 2 + N 3 )/(N 1 - N 3 ) . 


( 1 -b) 


The tacit assumption made of the N^ appearing in these 
expressions is that they are directly from the source (i.e., 
no scattering has contributed) and that the BATSE detectors 
are perfect. Any error source*in the will, of course, 
give rise to uncertainties, 69 and 6 <|> , in the source 
location. Thus, in utilizing counts to make the source 
location determinations, it is imperative to identify all 
possible sources of error in the recorded counts and to 
compensate for these in the location calculations. 

GRO will orbit at an altitude of some 400 km., well 
above the Earth's atmosphere. It might be anticipated that, 
in the course of a given burst, a significant number of 
gamma rays might scatter from the Earth's atmosphere and 
enter the Earth-facing BATSE detectors. Such contributions 
would obviously produce misleading data and corrections for 
the scattered contribution must be made. An extensive study 
of the atmospheric scattering problem has been previously 
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undertaken by Moris 3 using Monte Carlo calculations to 
examine probable photon interactions with the atmosphere. 

The results of that study were presented in graphical form 
in a way indicating the relationship of the ratio of 
observed-to-incident flux from a given zenith and satellite 
location as a function of gamma ray energy. 

The objective of the present work is to develop a model 
for the scattering contribution from which analytical 
expressions for the anticipated scatter contributions can be 
derived. The contributions would then be used to correct 
Equations (1) and, hopefully, lead to a refinement in the 
source location computation. In addition, analytical 
expressions offer the advantage of presenting statements 
which can be more readily interpreted than numerical 
results; it is thus hoped that this work can shed further 
light on the difficult scattering problem. 
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GAMMA RAY INTERACTIONS WITH THE ATMOSPHERE 


Gamma rays entering the Earth's atmosphere may be 
expected to interact with the electrons of that medium 
(bound primarily to nitrogen and oxygen, whose abundances 
are roughly 80% and 20%, respectively). The three principal 
interactions will be photoelectric absorption, Rayleigh 
scattering, and Compton scattering. The first of these will 
serve to attenuate the burst intensity as it propagates 
through the atmosphere, while the two types of scattering 
processes produce direction changes, along with an energy 
degradation in the case of Compton scattering. 


Suppose a gamma ray enters the atmosphere between the 
GRO orbit and the Earth's surface as depicted in Figure 3. 
Consider that portion of the incident flux, I_, which 
propagates through a column of cross sectional area A = 1 
cm^ and depth d, at which point scattering finally occurs. 
A typical photon energy in the range of interest in the 
BATSE may be taken as 100 keV. Relative to this energy 
electron binding energies are small 4 and the scattering 
electrons may be approximated as being free; thus the 
scattering may be regarded as Compton, with total cross 
section a * Since the total cross section represents the 
probability of a scattering event per electron, the photon 


certainty: no = l r meaning that 1/a electrons must be 

encounterd by e a c photon in the column or Figure 3 before 
scattering is certain to occur. If the atmosphere is 
regarded as being composed of two types of. atoms, a and b, 
with densities n (r) and n b (r) the number of electrons 
within the collum will be 


n e = Z a / n a (?) dv + Z b / n b (?) dv ' 


( 2 ) 


Z being the atomic number. Eguating this to 1/cr in essence 
gives a numerical value to the integrals since tne Compton 
cross section (as a function of energy) is well known. 

Next the photoelectric attenuation over this collum 
depth must be evaluated. The intensity diminishes by dl 
along some infinitesimal path length dg of the collum 
according to 
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dl = - y Id £ , 


(3) 


where y is the attenuation coefficient and is related to the 
electron population, Zn, and absorbtion cross section, a , 
in product form. With variable densities of two types ^ 
and a here representing the photoelectric cross section, 
the intensity diminishes as 


dl = - I 


tz a n a (f) + 2 b n b<*” 


a d£ 

y 


(4) 


Multiplication and division of the right hand side by A 
gives 


I = I Q exp {- cr Y [Z a / n a (?) dv + Z b / n b (?) dv] } . (5) 


The integral part of the argument is just n e , that 

is 1/a c , giving the attenuation of the incoming beam along 

the scattering path as 


I = I Q exp [- a y /a c ) 


( 6 ) 


If the scattered beam is imagined to undergo the same 
magnitude of attenuation in its (outward) propagation, the 
inclusion of a factor of 2 ion the argument of Eguation (6) 
will relate the attenuated intensity which might reach GRO 
in terms of the incident intensity undergoing a single 
scattering event. The specific form of the energy-dependent 
absorbtion cross section is given by 


a 

y 


4/2 

( 137 ) 4 « 7/2 


(7) 


where® Z here will be taken as an average for the 80% 
nitrogen (Z = 7) and 20% oxygen (Z = 8) composition, a is 
the Thomson cross section' (= 6.65 x 10" 25 cm 2 ) and a ?s the 
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ratio of photon energy to the electron rest energy. The g 
energy-dependent Compton cross section can be expressed as 


(8) 


where f(e) is a decreasing function of energy. Finally, 
putting together these results relates the (single) 
scattered intensity that might reach GRO and the incident 
intensity as 


I - I exp {- (6.2 x 10“ 4 )/f (e)ct 7/2 } . ( 9 ) 


Elementary considerations^ of a scattering situation in 
which a beam of intensity I Q is incident normally onto a 
thin target material anticipate that a detector some 
distance r from a scattering center will receive flux 
scattered between 0 and 0 + d0 (measured from the incident 
beam direction) given by 


I' 



, H rr . 

I r* n I 


( 10 ) 


where da/dg is the differential scattering cross section 
(per electron) for the particular scattering event. 
Scattering occurs symmetrically about the incident beam 
direction and it should be noted here that if the beam were 
oriented at some angle t, with respect to the scattering 
surface, the normal component would be I Q cos 5, and the 
cos £ factor would be included on the right hand side of 
Equation ( 10 ). 

Since the BATSE detectors count photon numbers rather 
than energies it is convenient to describe the scattering on 
terms of numbers. If a certain portion, dl^, of the inci- 
dent flux is due to dN photons at energy e Q , dI Q = dN Q e o . 
Similarly, if dN' is the portion of the total number 
scattering at energy e' , dl 1 = dN' e'. Equation ( 10 ) can 
then be expressed, allowing for non— normal incidence, as 


dN' 


dN Q cos n 




l dn J 


(ID 
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The two types of scattering of interest here are, as 
previously mentioned, are Rayleigh and Compton. Rayleigh 
scattering is elastic, while in the Compton case the ratio 
of energy of the incident to scattered photon energy can be 
simply expressed in terms of the scattering angle: 


e 


o 


1 + a(l - cos e) 


(12) 


Subscripting the 0-dependent differential scattering cross 
sections of the two types for purposes of distinction qives 
then 


dN' 


dN cos r 
o s 


2 
r . 


r r dq( 0 ) 

U d n 1 


R 



/• dq ( 0 ) 
1 dg 


) 1 


f 


(13) 


where the energy ratio is given by Equation (12). 

The detailed forms of the differential scattering cross 
sections depend, of course, upon the specifics of the 
interactive processes operative in the particular scattering 
event under consideration. The forms appearing in Equation 
(13) are given by lu 


and 


11 


r— i 


R 


-8 z3<J n 
(1.3 x 10 B ) — ^ 

a 


^1+COS 9^ 

. 3 ©J 
2sin j 


(14) 


r o 2 (l+cos 2 0 ) 

dfl j “ 3 {1 + 5 

C 2 [1+a ( l-cos0 ) ] (1+cos 0) Il+a(l-cos0)] 


(&) ■ 


a 2 (l-cos0) 2 


}. (15) 


Equation (15) is the Klein-Nishina formula and r n = (3cr /in 

1S 7 re ^ erre< ^ ^he "classical electron radius"* whose°value 

is ' 2.82 x 10~ 13 cm. 
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THE SINGLE SCATTER MODEL 


-.e-e — j — — * r\c* a 

C1LCLLO Wit LUC 


V> -» v*/^ v> >m r 

UUJ. Vi A 1. Ujr 


Atmospnerxc SUcill ksl my 

experiment have been discussed in some detail by Langer et 
al. Their analysis indicated that, given the atmospheric 
composition* , 30 keV photons should scatter at a height of 
about 40 km above the Earth's surface, while 500 keV photons 
scatter near 30 km. The enerqy range of interest to BATSE 
(maximum detector efficiency range) is, roughly, 25 keV to 
250 keV. The above scatter altitudes are thus relevant to 
the present study and their slight difference in comparison 
with the 400 km GRO orbit suggests representing the 
scattering as taking place from a sphere of radius ~40 km 
larger than the Earth's radius. 


Limiting the present study to single scatter 
contributions, all of which are imagined to come from a 
spherical surface of radius about 6411 km, means that any 
photons scattering into the BATSE detectors must have 
scattered from a point within the spherical cap depicted in 
Figure 4. The scattering region (cap) is defined by imaging 
a line extending from GRO tangent to the scattering 
sphere. Allowing that tangent to sweep around the vertical 
from GRO to the sphere's surface will define a cone whose 
intersection with the scattering sphere gives the region 
from which single scatter events might effect BATSE data. 
That GRO is at altitude 360 km above the scattering region 
means that the cap might be reasonably approximated as a 
flat disc. 


Consider then the situation shown in. Figure 5-a. N 
photons (per cm^ per sec) are imagined incident onto a plane 
at some angle ?. The spatially uniform flux over a disc of 
area Ap will give rise to dN cos £ striking the area of the 
washer of radii r' and r' + dr'; dN Q can be expressed in 
terms of N Q as 


dN 

o 



2irr ' dr ' 



1 • 


(16) 


All points within the washer are equidistant, r, from GRO 
which remains at altitude h above the plane. As previously 
mentioned scattering will be axially symmetric about the 
incident axis. Thus, only that small fraction of photons 
scattered through the required 0 and into A<j>/2ir can 
contribute to detector counts. But A<fr is the axial angle 
subtended by GRO dimensions, A<f> = Aw/r, and this must 
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multiply, the dN' in order to obtain the possible number of 
photons scattered into the GRO direction. 

The above results are now to be assembled in order to 
compute the anticipated numbers of photons scattered into 
the BATSE detectors. The number possibly scattered into GRO 
is 


dN' 


= AW_ 
GRO 2 tt r 


dN' 


(17) 


dN' is that number scattered (per electron) and is computed 
by using Equation (16) in Equation (13) with the attenuation 
factor given in Equation (6) multiplying (i.e., diminishing) 
the incident number of photons per cm 2 per sec. These steps 
give 


dN GRO 


= N cos r 
o 


AW 


e -'V'c> f 


r'dr' 



That the differential scattering cross sections refer to 
scattering per electron means that <3N' GR0 must finally be 
multiplied by the total number of scattering electrons 
imagined to reside on the scattering disc. Assuming the 
uniform distribution 6 cm - over A^ means that <$A D 
scattering electrons reside on the disc. The total number 
of photons that can scatter into GRO is the product 
of 6 A d and Equation (18): 


dN GRO 

TOT 


N cos 
o 


?5AWe 


- ( o/° c > 


[( 


da 

dn 


) ♦ 
R 


c£) 


(§§> 1 


r'dr' 


(19) 


The relevant geometry of Figure 5-a has been 
highlighted in Figure 5-b. It is evident that 


r 1 

tan 4 > = ; 


( 20 ) 


further, since 
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0 + <|>+5 = nf 


( 21 ) 


it follows that 


*1 

h 


tan ( e + e ) . 


( 22 ) 


In order that r'/h > 0 it is apparent that -|- < ( 0 + c ) < ir . 
Equation (21) is equivalent to 


0 “ IT “ £ “ (f> 


0 __ corresponds to the minimum value of <(> while 0 _._ occurs 
wh'in <|> takes on its maximum value. The minimum vai‘6e 
of <|> is obviously zero, while the maximum occurs when r' 
takes on its maximum value. The latter is simply the radius 
of the scattering disc which examination of Fiqure 4 shows 
to be sin (18.5°) , R A being the 6411 km radius of the 
scattering sphere. Thus, 


_ x R A sin(18.5°) 

0 = tt - c ; 0 • = w - r - tan f r 1 

max ^ min * L h 1 


(23) 


It is a simple matter now to convert the integration 
variable of Equation (19) to the scattering angle. It 
follows immediately from Equation (22), along with the 
Pythagorean relationship of r and r' that 


r'dr ' 


= K sin 


(0 + o d0 . 


(24) 


Substitution of Equation (24) into (19) gives 
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dN GRO 


N q 6AW "(cr /a ) 


Y c {cos 2 ?[(^) Sinede + (^t)(§§) sinede] 


(25) 


+ coscsin5 U 35 ) cosede + coe9<J «H 


Evaluation of the four integrals arising from Eguation 
(25) is necessary to now complete the assessment of the 
number of photons that might be expected to scatter into the 
BATSE detectors. The integrals are: 


/ (3§) sinede 

e . u r 

min 


3 e 

q a Z max . . 2 n 

(1.3 x 10 ) ■ °- f — — f - sinede 

a e . 2sin J 4 
min 2 


max . 

f cosede 


0 min 


2 9 

_o a Z max., 2 n 

(1.3 x 10 ) — 5 — f = — — cosede (27) 

<* e m . 2sin' 3 | 
min 2 


max e , 

/ sinede 

e min 


l+cos e 


O , , 1 + CO 

16 " LJ U+«<1- 


[1+a ( 1-cose) ] 


+ -gl * - 1 — c - ° - qa > 2 [ cosede 

[ 1 + a ( 1-cose )] J 


(28) 


9 max e_ . 3 a 6 nax . 2 

/ (7?) (§§) cosede - jg ! ( ltcos 9 2 

9 min c 16 ” 9 mln 1+ad-cose)] 2 


+ e 2 ( 1-cose) 2 , oos6d9 

[ 1 + 0 ( 1-cose ) 


(29) 
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The integrals appearing in Equations (26), (27), and (28) 
are easily reduced to algebraic forms, while that in 
Equation (29) is not so simple. The results are the 
following: 


0 „ 3 

max , _p a Z . 

/ (SJ) sinede = >/2 ( 1 . 3x10 S )-^- f 

0 . R 

mm 


a 3 \| 1-cose 
0 . 


(30) 


+4 J l-cose - 1-cos 0) 3 ^ 2 ] 


max 


0 . 
mm 


max^ /T/i •> m-8 s CT o Z f 7 r 4 1 — cos0- 2-| 

/ (do) D COS0d0 = (1-3x10 ) in [ : ^ TT __ rT ] 

0 . R a 2<J2 

mi n 


(31) 


0 

-JT+cos? [(4+|cos0) + 1 _c'ose 1 > 


max 


6 . 
mm 


f maX f^)(f§) ^ (l-cos») - ,,! 2 ° Wl) 

0 . e u c 
mm 


1 6 it 1 2 

a 


a [l+o( 1-COS0 ) ] 

(32) 


2(l-cos0) + 


[ 1+a ( 1-cos 0 ) ] 


3_ 0 . 

- K ( l+2a ) in [l+a( 1-COS0 ) ] } 

] 1 a 0 


max 


mm 
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0 


max e 


min 


(T?Ha§) 


3 a . 

T57 W (2+2a 


ct ) 0 + — ■ j sine 


a(l+q) ( 1+COS9 )sine 
2 ( l+2a) 2 [ l+<x( l-cose ) ] 2 


l ( lla^+llq^-30a^-48a 2 -24q-4 )sin 9 
2 ci 3 (l + 2 a) 2 [l+a( l-cose ) ] 


(33) 


0. 


+ tan -1 (^T+ 2 cT tan A)} 


max 


min 


Evaluation of these results over the limits given in 
Equation (23) gives, upon multiplication by the factors 
indicated in Equation (25), the inteqrated form of Eauation 
(25) which represents the anticipated number of photons 
scattered into the detector direction. It is apparent that 
the Rayleiqh contribution is wholly ignorable. 

The angle ? will not be known in general and so must be 
expressed in terms of more appropriate quantities, namely 
satellite location angles, (e , a), and burst If cation 
coordinates, ( 0 , <j> ) , computed from Equation (1). 

Reference to Figure 6 shows that ? is defined locally as the 
angle between the burst direction and the position vector of 
GRO. Represent, for the moment, the burst location in terms 
of geocentric angular coordinates (0 , <|> ). The incoming 
burst flux can then be represented a£ ^ 



Let be the GRO location at this time: 

*o = R 0 (sine o cos<|> 0 i + sine^in^ j + cos 0 Q k). (35) 

The angle 5 can easily be seen as 
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K -*n 

COS ? = (r”)*( 7) “ COS 5 (® Y ' ♦y' 0 O ' ^ 0 * 

o 

The burst location, however* is*inade not in terms of 
(0 , <j, ) but in terms of (0 , * ). Use of the known 
orientation of the satellite axes relative to the geocentric 
axes enables 0 and 4 to*be expressed in terms of 
calculable functions Xf 0 and <|> .* Us$ of the latter 
transformations give cos? = cos?(0 , 4> 7 9 Q , <t> ) and 
expressing ? in terms of these coordinates in the integrated 
form of Eguation (25) finally gives the scattered 
contribution to GRO in terms of burst location and satellite 
coordinates . 
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DISCUSSION 


The objective of the present work was to develop 
analytical expressions for estimating the atmospheric scat- 
tering contribution to detected photon counts in the course 
of a burst event. Having such expressions in terms of source 
and satellite location coordinates would, presumably, furnish 
a more descriptive picture of the scattering effects than 
would be offered by purely numerical results. Although the 
objective seems to have been achieved in principle, by means 
of a model calculation, the cumbersome nature of the final 
expressions does not allow a guick and easy grasp of the pic- 
ture to be made. Their most important contribution might be 
made by their use in examining various important cases not 
covered in the Morris' 3 study. To this date only limited num- 
erical evaluation of the results have been carried out. In 
particular the shape of the back-scatter curve as a function 
of photon energy for c = 0° and an infinite disc (an approxi- 
mation which enormously simplifies evaluation over the 
limits) has been investigated and found to be in good agree- 
ment with the Morris-type curves. Inputting the final forms 
to a computing facility should produce curves for any desired 
situation . 

The flat disc model used here was, obviously, for ease 
in calculation. It indicated that such calculations are 
indeed possible and so a redevelopment of the results for a 
spherical cap should be managable. The fact that, at a 400 
km orbital altitude, GRO will see only about 2.75% of the 
Earth's surface seems to make the disc model reasonable, but 
not perfect. In particular, a spherical cap could offer a 
"shadow region" to a burst, as the curvature would hide that 
part of the scattering region away from the burst direction; 
the flat disc, of course, offers no such feature and thus 
might suggest some additional contributions to scattering 
which would not actually occur. 

The model discussed also considers only single scatter 
contributions. The Rayleigh contribution, as already noted, 
is completely negligible for photon energies of interest here 
in the case of interaction with the relatively light elements 
nitrogen and oxygen. Only Compton scattering would thus 
significantly contribute any photons. Heitler and Nordheim 14 
are credited as being the first to consider double Compton 
scattering and found the double cross section to be down from 
the single cross section by a factor of 1/137. The error 
made by ignoring double scattering would thus seem to be 
acceptable. The fact that single scatter events dominate 
would seem to offer validity to the nation introducted here 
that the only significant scattering contributions are those 
coming from a cap region immediately below GRO. 
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Figure 3. Photon Propagation Through a Column 
of Atmosphere Below GRO 



Figure 4. Single Scatter Region 
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Figure 6. Schematic of Position Relationships 
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